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Abstract 

Qs ', We study the properties of the spin dependent one body density in momentum 

space for odd-A polarized deformed nuclei within the mean field approximation. We 
derive analytic expressions connecting intrinsic and laboratory momentum distribu- 
tions. The latter are related to observable transition densities in p-space that can 
^ I be probed in one nucleon knock-out reactions from polarized targets. It is shown 

C ■ that most of the information contained in the intrinsic spin dependent momentum 

distribution is lost when the nucleus is not polarized. Results are presented and dis- 
ked I cussed for two prolate nuclei, ^^Ne and ^^Mg, and for one oblate nucleus, ^^Ar. The 
V^ ' effects of deformation are highlighted by comparison to the case of odd-A nuclei in 
the spherical model. 
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I. Introduction 

Momentum distributions in nuclei are a subject of increasing interest at both nucleon and 
parton levels. Reviews of theory and models on nucleon momentum distributions can 
be found in Refs. |]l], 0. Experimental information on nucleon momentum distributions 
can be gained from one nucleon knock-out in nuclear photo-absortion, as well as by 
A(x, x'N)B reactions'^ -where a hadronic or leptonic beam, x, scatters and knocks out 
one nucleon from a nuclear target, A. As pointed out in Ref. [^ the interpretation of 
information on momentum distributions is not free from ambiguities because, at present, 
neither experimental data nor theoretical considerations are sufficiently comprehensive. 
Yet fairly reliable information is already available on nucleon momentum distributions 
from electron scattering in the quasi-elastic regiono, as well as on parton distributions 
from lepton scattering in the deep inelastic regionB. Most of the work on complex nuclei 
has been devoted to spherical nuclei, but in the last years several questions have emerged 
concerning the role of nuclear deformation in quasi-elastic electron scattering and in 
deep inelastic scattering. In the latter, it has been suggested to use polarized deformed 
nuclei to study the degree of transparency in different directionsO ; in the first it has been 
questioned whether coincidence measurements may be sensitive to nuclear deformation. 

Investigations of quasi-elastic electron scattering from even-even deformed nuclei have 
been carried out both theoretically and experimentally. Starting from mean field calcu- 
lations with increasing degree of sophistication, we studied the properties of momentum 
distributions in the intrinsic ground state of several even-even axially symmetric deformed 
nuclei. We investigated not only effects of deformation but also effects of pairing and of 
short range correlationsB & a . In a search for deformation effects, coincidence Cfije'p) ex- 
periments on spherical and deformed Nd isotopes were performed at NIKHEFli3. More 
recently we have also initiated studies of deformation effects in quasi-elastic electron 
scattering from polarized deformed nuclei EJ. With progress made in developing polar- 
ized beams and targets, as well as in polarimeters designs (in addition to that on CW 
linac facilities), the study of spin degrees of freedom is becoming a major theme in Nuclear 
Physics either at low, intermediate or high energies. In this paper we attempt to lay a 
basis for such studies in odd-A deformed nuclei. 

In general the one-body density (either in r-space or in p-space) is a two by two matrix 
in spin space. For closed shell (spin saturated) nuclei the one-body density is proportional 
to the unit matrix in spin space, i.e., the momentum distribution is independent of spin. 
The same is true for even-even axially symmetric deformed nuclei in the mean field 
approximation due to time reversal invariance of the intrinsic ground state. However, for 
an odd-A deformed nucleus the intrinsic ground state is no longer time reversal invariant 
and the one-body density matrix is in general non diagonal in spin space. Hence scalar 
and vector momentum distributions can be defined in the intrinsic frame, that contain 
all the information on the intrinsic spin dependent momentum distribution. This paper 
is devoted to the study of these new momentum distributions. 

It is important to realize that in order to access to this new information one needs 
to relate the observable spin dependent momentum distributions of the polarized nucleus 
in the laboratory frame with the intrinsic scalar and vector momentum distributions. A 
central issue of this paper is to establish such relations. 



We study properties of scalar and vector momentum distributions within the self- 
consistent mean field approximations assuming axially symmetric and time-reversal in- 
variant mean fields. To this end we present and discuss results of density dependent 
Hartree-Fock (DDHF) calculations for two nuclei ^^Ne and '^^Ar that have ground states 
with equal spin and parity (3/2+) but have quite different intrinsic structure. Results for 
^^Mg are also discussed. The self-consistent mean fields for ^^Ne and ^^Mg are prolate 
whereas that for '^^Ar is oblate, and the i^-values are different for each nucleus. The role 
of deformation is further highlighted by comparison to the limit of a spherical mean field, 
i.e., to the case of odd-A nuclei with a single-nucleon outside closed shells. 

The paper is organized as follows: in Section II the intrinsic scalar and vector momen- 
tum distributions are defined and their properties are studied. In Section III we define 
scalar and vector momentum distributions for the polarized nucleus in the laboratory 
frame, we establish their connection with the intrinsic ones, as well as with observable 
transition densities in p-space. In this section the spherical limit is also discussed and 
comparisons with results on the deformed nuclei are presented. Section IV summarizes 
our conclusions. 

II. Intrinsic Momentum Distributions 

To study spin dependent momentum distributions in axially symmetric deformed nuclei 
our starting point are the single-particle Hartree-Fock (HF) solutions in coordinate space, 
which are characterized by the angular momentum projection along the symmetry axis 
(Qj) and by the parity (ttj). The HF wave functions are expressed in the intrinsic system 
with the z axis the symmetry axis and the (x, y) plane the symmetry plane. The states 
i, i degenerate in energy are written asliS 

^Ar.s^q) = xM [<l>^\r^,z)e'^^-x+{s) + <f>t\r±,z)e'^^\4s)] , (1) 

Hr,s,q) = xM [-4'\r^,z)e-'^''\+is) + 4'-\r^,z)e-'^^-X-{s)] (2) 

where Xgiil)^ X±{^) are isospin and spin functions, i is the time reverse of i, A"^ = 
r2j±l/2 > 0, and r±, z, tp are the cylindrical coordinates of r. The wave functions (f)\ are 
expanded into eigenfunctions of an axially deformed Harmonic Oscillator, characterized 
by the set of quantum numbers {a} = {ra_L, ra^. A, S}, with expansion coefficients C^. 
Taking the Fourier transforms of Eqs. (1) and (2) we get the single-particle HF wave 
functions in momentum space 

Up,s) = 4>\''\p±,Pz)e'''^'''x+is) + 4>['\p±,P.)e'''^^\4s), (3) 

$j(p, s) = [-4>l-\p^,p,)e-'^^^\+is) + 4>\''\p±,Pz)e-''''^'xAs)] ^^ (4) 

with p±, Pz, ^Pp the cylindrical coordinates of p. Here and in what follows we omit isospin 
to abbreviate the notation. The wave functions (j)} are given by 



^(T) 



(pi iP±,Pz) = J2 Cl,±(i>a±{p±,Pz 



(5) 



n±,«z 



where a^ stands for {n±, Uz, A = A^} and the sum extends over basis states with major 
quantum numbers (A^ = 2n± + nz + A) even or odd for ttj = +1 and ttj = —1, respectively. 
In the cyhndrical basis 



<t^a±{P±,Pz 



-t)'^^/rci .± 
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(6) 



(7) 



(8) 



with L„^ and i/„^ associated Laguerre and Hermite polynomials respectively, rj = p\/(3]_, 
^ = Pz/Pz and f3±, j3z the inverse of the harmonic oscillator length parameters {f3± = 
{mu±/hY^'^, Pz = {mcUz/hy^'^). For later use it is convenient to expand the momentum 
dependent single-particle wave functions in a spherical basis. We write 



Hp,s)=j:yi^'ip,s)<i>np) 



(9) 



l>j(p, s) = ^(-l)^-^-3^,/'(p, s)0f (p) 



(10) 



with 



and 



y^s) = T.^^A- S|j m)Yt{p)x^{s) 



AS 






:iii 



(12) 



Rl{p) = I dQpY/'*{Q,)MP±,Pz)e'''^'' 
(2£+l)(£-A)!T'/' 



'-^)''^'\/7r,, 



2{i + A)\ 



sin 6pd6pP^ {cos Op)ip'^^ {p sin Op)ilJn^ [p cos Op) . 

(13) 



We represent by \E'fc the intrinsic ground state of an odd-A nucleus with odd nucleon 
in state k. $fc represents the wave function in momentum space for the state k, and tt^, 
Vtk = K are its quantum numbers. The intrinsic one-body spin dependent momentum 
distribution is then given by 



Mss'(p) = (^fc|apE«ps'l*fe) 
= T^^i [(Xs'|$i(p, s)){^i{p, s)|xe) + (Xs'l^ib, s))($i(p, s)|xs) 

+ (xs'|$fc(p,s))(l>fc(p,s)|xs) (14) 

with vf occupation probabilities. The last term in this equation gives the contribution 
from the odd nucleon, while the first term gives the contribution from the nucleons in the 
even-even core. Decomposing this matrix into a scalar (M) and a vector (M) in spin 
space 

M = Tr[M{p)] (15) 

M = Tr[M{p)(T] (16) 

it is easy to show that the even-even core contributes only to the scalar momentum 
distribution, 

M(p) = Y.^v^ {\4>\''\P^,P^)\' + \4>t\p±,P^)\') 

+ \~^^:\p^.p.)?+\~A'\p±.p^)\\ (17) 

while the vector momentum distribution is made up entirely by the odd nucleon, 

MM = \4>i^\p±,P.)\' - \4>i~\p±,Pz)\' (18) 

M,(p) = cos(/?pMx(p±,P.) (19) 

My{p) = simppM^^^p,) (20) 

with 

MM±,P.) = 24>i^\p^,pM'\p±.Pz) . (21) 

In previous worko El we studied momentum distributions in even-even deformed nu- 
clei. For spin-saturated even-even nuclei the vector momentum distribution is zero in the 
intrinsic ground state and only the scalar momentum distribution remains. To connect 
with the notation in Refs. [^, ^ we denote by n{p) the contribution from the even-even 
core to M{p), and by nx{p) its multipoles. The contribution from the odd-nucleon is 
denoted by n°'^'^{p) 

M{p) = n{p)+n'""'{p) = Y. ^a(cos^p) [nx{p) + nf'ip)] = ^ P^icos dp)M x{p) . 

A=even A=even 

(22) 
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A deformation parameter in momentum space (3'^ was introduced in Refs. p, ^ defined 
in analogy to the standardEJ definition of the mass quadrupole deformation parameter 

V 5 Jn{p)p'^dp 

For odd-A nuclei the same definition applies replacing n{p) by M{p). 

In what follows we show scalar and vector intrinsic momentum distributions of ^^Ne 
and ^'^Ar. The results presented have been obtained using the SKA effective interactionEil 
and the McMaster version of the deformed Hartree-Fock codeO that follows closely the 



procedure of Ref . |]T2| . The HF results for binding energies as well as for the charge and 
mass quadrupole moments and r.m.s. radii are summarized in table I. Also given in the 
table are the results for the even-even isotopes ^°Ne and ^^Ar. Important differences 
between the structures of Ne and Ar isotopes are that the former is prolate while the 
latter is oblate and that the unpaired neutrons in ^^Ne and in ^^Ar sit in orbitals with 
Qk = 3/2 and with fi^ = 1/2, respectively, while both nuclei have J = 3/2 in their ground 
states. Also given in table I are the moments of inertia and the decoupling parameter 
(defined as in Ref. |1^). It can be checked that with these parameters the angular 
momentum sequences of the ground state bands in ^^Ne and in ^^Ar are well describedO 
when compared to experiment B. The moments of inertia have been calculated using the 
cranking formula. 

In Figure 1 we show in three-dimensional plots the contributions to the scalar mo- 
mentum distribution in ^^Ne (top) and ^^Ar (bottom) from the core (left) and from the 
odd nucleon (right). As already mentioned, the odd nucleons in ^^Ne and ^^Ar have 
Qk = 3/2 and 1/2, respectively. Clearly the odd nucleon contributions are quite different 
in ^^Ne and in ^^Ar, this is due to the different K values and mean field deformations. 
We shall come back to this point later on. It is also interesting to remark the difference 
between odd nucleon and core nucleons contributions. An important difference between 
the momentum distribution of the core {n{p)) and the momentum distribution of the odd 
nucleon {n°'^'^{p)) is that in the first case the monopole term, rio(p), is dominant while 
in the second the A = 2 and higher multipoles are comparable to the monopole term. 
In Figure 2 we show the A = 0, 2 multipoles of the even-even core of Ne and Ar. It 
is worth pointing out that the isotropy condition found in previous work on even-even 
deformed nucleii' B is also met here for the even-even cores, in the sense that when we 
compute the deformation parameter (3^ we get /3^ ^ (see table II), whereas the corre- 
sponding parameter in r-space is large, consistently with the experimental deformation. 
Experimentally^' Il2l the quadrupole moments are 51.4 ± 3.8, 58.5 ± 2.6, and —38.5 ± 21 
fm^, whichcorrespondto/?'" values of 0.45 ±0.03, 0.51 ±0.02, and -0.15 ±0.08, for ^^Ne, 
^°Ne, and ^^Ar, respectively. 

Actually, as seen in table II, (3^ is close to zero not only for the even-even, but also 
for the odd-A nuclei. This important property follows from the major shell admixtures 
in the single-particle wave functions that result from the self-consistent mean field® B. 
Interferences between different major shells have opposite sign in r-space and in p-space, 
and make it possible to reach configurations with large equilibrium deformations and 



minimum kinetic energies. To illustrate the role of N-admixtures in momentum space we 
show on the right hand side of Fig. 2 the results obtained with the Nilsson model when 
the Nilsson hamiltonian@ is diagonahzed in the entire basis space (N< 7) and when it 
is diagonahzed within each major shell subspace (the x-, l^ parameters have been chosen 



as in Ref. [^). Clearly, the results for n2{p) are substantially different in both cases. 
When AN admixtures are taken into account the results are similar to the ones obtained 
in DDHF: the oscillations in n2{p) are such as to minimize the (3^ value. When AN 
admixtures are neglected the A = 2 multipole of the density in momentum space has the 
same shape as the A = 2 multipole of the density in r-space, and the (3^ value turns out to 
be the same as the deformation parameter in coordinate space j3'^ (see table II). A similar 
situation is found for the A = 2 multipoles of the odd nucleon momentum distributions 

The vector momentum distribution components Mz{p) and Mj_[p) are also shown on 
3-dimensional plots in Figure 3, where the results on the top correspond to ^^Ne and the 
results underneath correspond to ^^Ar. As seen in this figure, the momentum distributions 
in ^"^Ar have a richer structure than those in ^^Ne. This particularly rich intrinsic structure 
is characteristic of deformed odd-A nuclei with K = ^. A characteristic of nuclei with 
K > ^ is the nodal line along the axis p± = (see the top plots of Fig. 3). As a rule, the 
structure tends to become simpler for the higher K values. As an example, we show in 
Figure 4 the momentum distribution of ^^Mg, that has K = ^. In this case {Qk = K = j), 
Mz{p) = 71°^*^ (p) and M±{p) = 0. The function ra°'^'^(p), represented in the top of Fig. 4, 
shows a simple structure corresponding to a 99% d^/2 wave function. At the bottom of 
Fig. 4, one can also see the core momentum distribution in ^^Mg which is very similar to 
that of 21 Ne. 

To understand the structure of these figures it is convenient to express the densities 
in momentum space in terms of the spherical components of the single-particle wave 
functions. Using Eqs. (|^)- ([T2|) one can show that 



^i + kk ii^-^'l 






m4p^,p.) = ^ E E alp4cos^,)(J ^ o)n°,,(p) (25) 

L=even A=odd \ / 




47r .f:-„_^.. /rr^iu;^'^"""^ni -1 0; 



L=evenA=odd y {L + 1)L 

where we have introduced the momentum dependent density 

nUp)=EE0gF(p) ^' L i \i-ir^''^''{jKf -K\\0) (27) 
ij t-'i' [1/2 1 I/2J 
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with 



QlT(p) = ^npW(py^'Ln'r j^ -] (2j 



We use the notation a = y/2a + 1. Eq. ( pSD also apphes to 9/^ ■' (p) with k replaced 

by i and L by A. (j)i^{p) {4>i{p)) are the spherical components of the odd nucleon (core 
nucleons) wave function defined in Eqs. (0) and (|l^). Note that under interchange of p^ 
by —pzi i-e., reflection through the plane perpendicular to the symmetry axis, M{p) and 
Mz{p) are unchanged, whereas M±{p) changes sign. 

The fact that each individual single-particle wave function contains in general more 
than one ij component, makes it possible to have contributions to Eqs. (p^)-(|27|) coming 
from several £j-components, as well as contributions from interferences between different 
ij waves. This is an important feature that would not be possible in the spherical limit, 
i.e., within the extreme spherical shell model. The lower the K value is, the lower ij- 
values are allowed in the single-particle wave functions, producing a richer structure in 
the momentum distributions (compare the results for ^^Mg, ^'^Ar and ^^Ne in Figs. 1-4). 
In these sd-shell nuclei the dominant ij components of the outermost nucleons are i = 
and £ = 2, though higher components are also present in our single-particle HF wave 
functions. In particular, the i = wave is allowed for the K = 1/2 wave function while 
it is not for K = 3/2, 5/2, giving rise to the richer structure observed in ^^Ar. 

III. Momentum Distribution Spin Matrix in Labora- 
tory Frame. Polarization Considerations 

To study spin dependent momentum distributions in laboratory frame we need to consider 
polarized nuclei. We consider a nucleus in its ground state with angular momentum J 
completely polarized in a given direction P* defined by the angles Q* = {6*,(f)*) in the 
laboratory frame. The two by two matrix 

M!;i}{p,n*) = {JJ{n*)\a^,ap,,\JJ{n*)) (29) 

measures the probability to find a nucleon with momentum p and spin projections a, a' 
in the laboratory frame. Our task here is to relate momentum distributions that may be 
observed in the laboratory to the intrinsic momentum distributions studied in the previous 
section. For this purpose it is also useful to consider partial contributions to Eq. ( [29| ) in 
addition to the total momentum distribution above defined. We note that Eq. (|29|) can 
also be written as 

Mi'Xp^n*) = Y.{JJmHJR){R\ap.'\JJ{^l) (30) 

R 

where the sum is carried out over a complete set of states R of the residual system with 
A-1 nucleons. In coincidence reactions of the type (x, x'N), where x represents a leptonic 
or a hadronic projectile and N a knock-out nucleon detected in coincidence with the 



scattered particle x', one can select transitions to discrete states in the residual system 
and, therefore, one can measure partial momentum distributions of the form 

Mj;ir'''\p,n*) = Y.{JJm\a^JJnMn){jRMR\apMJ{n*)) ■ (31) 

Mr 

In what follows we study the connection between the momentum distributions defined 
by Eqs. (ES) and (RTl) and the intrinsic momentum distributions. 



III.l Total Momentum Distributions 

To compute the total spin dependent momentum distribution defined in Eq. (|29|) we first 
note that the magnetic substate \JJ{Q*)) must be expressed in terms of states referred 
to the same quantization axis as the spin components a, a'. We take this quantization 
axis to be the z axis of the laboratory fixed frame, 

\JJiQ*))=T.1^iiji^*)\JM) (32) 

M 

The convention for rotation operators and T> matrices is as in Refs. p^, |23[| . 

On the other hand, we have to write the angular momentum eigenstates \JM) in terms 
of the intrinsic ground state \E^fc of the odd-A deformed nucleus. For that purpose we use 
the Bohr-Mottelson factorization approximation a and we get after some algebra 






X 



+ 



^K -K o)(-l)'''''^io(^)[(*^l4E«ps'l*fc) + (-l)'(*^l4EapE'|^fc) 



,K K -2KI 



^M2x(^)(*fcl4s«PS'|^fc) + (-l)''P;_2x(^)(^fcl4s«PS'|*fc) 



(33) 



The same result is obtained in Projected HF to lowest order in (J^)~^ (see Refs. |T^, ^). 
We note that the first two terms within the brackets are the intrinsic spin dependent 
momentum distributions Mss'(p) ^^^ {—l)^~'^M^j._j.,{p), introduced in Eq. ([T^) (see 
also Eqs. (^) and (^)), whereas the two last terms are new. The latter depend only on 
the intrinsic momentum distributions Mz{pz,p±) and M±{pz,p±), which in turn depend 
only on the odd nucleon wave functions (see Eqs. (|1^) and(^)). We define M (p, fi*) 
and MJ^'^\p,n*) (with a = 0, ±1) as the trace and the spherical vector components, 
respectively, of the momentum distribution spin matrix in laboratory frame for the fully 
polarized nuclear ground state 

M^'^^ (p, n*) = Tr (m^-^^ (p, fi*)) (34) 

M^) (p, n*)=Tr (m(^) (p, n*)a^ . (35) 
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It is a simple matter to show that the scalar momentum distribution M [p^Vt*) is 
proportional to the intrinsic scalar momentum distribution and that the a components of 
the vector momentum distribution depend only on the odd nucleon. Hence, also in the 
laboratory frame, the core nucleons contribute only to the trace M (p, fi*) and not to 
M (p, fi*). Taking the trace in Eq. (^) and using Eq. (p2D it is easy to show that 

M^'\p,n*)= Y. Mx{p)Px{cose;)G{\;JK) (36) 

A=even 

where M\{p) are the multipoles of the intrinsic momentum distribution (see Eq. (0)), 
G(A; JK) is a geometrical coefficient 

G(A;J/0^(-l)--i^P !j ^)(,^ ^^ ^). (37) 

and 9* is the relative angle between the momentum p and the polarization direction P* . 
The expression for M {p, Q*) is somewhat more involved, as one gets contributions 



from all the terms in Eq. (|33|) that correspond to different angular momentum projections 
of the intrinsic odd nucleon densities M{p). After a lengthy but straightforward algebra, 
using Eqs. (^ and (||), we get 

Mi'\p,n*) = VQ y: E E(-i)%'m>^.%(^p) fi , f: , _0 ^i5(p) (38) 

L=cvcnA=odd M \a II a ^ij 

with 

£i5(p) = G{\- JK) |n°,,(p) + i-^y^''^^^^^}^^^ZiP)] (39) 

and a = 0, ±1. The momentum dependent densities are given by 

KMp) = Y.Y. 0gf ^'(p) \^' L i\ {-lY'{j K 3' K\\ 2K) (40) 

^] i'f [l/2 1 I/2J 

and n° j^(p) as given in Eq. (|27D. Note that the signature dependent term, Il1\{p), only 
contributes for A-values satisfying 2K < A < 2 J, i.e., A = 2J when K = J. 

If we fix the polarization direction as the z-axis in laboratory frame, Eq. (^8|) reduces 
to a simpler angular dependence involving the direction of p only: 

Mi'\p,^* = 0) = 76 E E {-ir^Y^\n;) (l ^^ ^) ci%{p) . (4i) 

L=even A=odd v ^^^ \ J 

We then define M\ , M]"') and M^^ as the vector momentum distribution components 
along the polarization axis (/) and in the transverse directions (s, n). Denoting by (6'* </)*) 
the direction of p in laboratory frame with the z-axis parallel to P* we have 



M/^)(p,fi* = o) = ^ E E ALP,(cos^;) (J \ ^)/:i» 

L=even A=odd \ / 



(42) 
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Mi'^ (p, n* = 0) = -| (Mf - m'l'}) ^ cos ^;Mi^) (p, n* = 0) (43) 

My)(P,fi* = 0) = ^(Mf^^+Mii^) =sm^;M;^)(p,fi* = 0) (44) 
with the transverse component defined as 

Mi^)(p,l]* = 0) = -^ E E 7=^^^2(cos^;)(| i^^ o)^i». (45) 



L=cvcn 



A=odd \J{L + 1)L 



The momentum distribution in the laboratory frame for the unpolarized nucleus can be 
easily obtained by integrating Eqs. (|36|) and ( p8|) over the polarization direction. Clearly 
for an unpolarized nucleus only the A = multipole of the scalar momentum distribution 
remains: 

/ dn*Mi^\p, ^*) = , / dn*M'^^\p, n*) = ATiMx=o{p) (46) 



For A = 0, Eq. (|2|) reduces to 

M,=o(p) = n.=o(p) + nf^M = ^ E 2 [^.' + kk (^ - vf)] ^ |0f (p)p . (47) 

Hence in the unpolarized case one loses not only all of the information contained in 
the vector momentum distribution but also important information contained in the scalar 
momentum distribution. One misses the A > 2 multipoles of the core and odd nucleon 
momentum densities that, as discussed in Section II, carry important information on shell 
admixtures and on details of the internal dynamics which are crucial to attain equilibrium 
in open shell nuclei. 

Now we can easily compare the scalar momentum distribution in laboratory and intrin- 



sic frames. Comparison of Eq. (p2|) to Eq. (^) shows that the dependence on the direction 
of p relative to the intrinsic symmetry axis, is replaced in the laboratory frame by a sim- 
ilar dependence on the direction of p relative to the polarization direction, weighted by a 
geometrical coefficient G{\\ JK). This geometrical coefficient takes the value 1 for A = 
and decreases as A increases {G{\\ JK) = for A > 2 J). For ^^Ar and ^^Ne the only 
non-zero coefficient with A > is G{\ = 2) = 1/5, —1/5, respectively. 

For the vector momentum distribution we also see that the expressions for longitudi- 
nal and transverse components in the laboratory frame involve the intrinsic multipoles 
n° x{p)- These expressions are similar to the ones in the intrinsic frame (compare Eqs. (|3^) 
and (|45|) to Eqs. ( pSj) and (pB])), replacing the direction of p relative to the internal sym- 
metry axis by the direction of p relative to the polarization direction, and weighting the 
intrinsic A multipoles with the geometrical factor G(A; JK). However in the case of the 
vector components one has additional contributions for the A > 2K multipoles coming 
from the signature dependent terms (see Eqs. (^) and (^OD). 
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III. 2 Partial Momentum Distributions 

In this section we consider partial contributions, as defined in Eq. (0), to the spin de- 
pendent momentum distribution previously discussed. These partial contributions play 
a central role in coincidence A(x, x'N)R experiments where transitions to discrete states 
Jr of the residual nucleus are selected. In these processes the dependence on the nuclear 
structure of the cross section contains as basic ingredients the spin dependent transition 
densities in momentum space M^^T defined in Eq. ( pi]) (see Refs. |]ll], ^). 



As seen in previous sections, where we considered total momentum distributions, only 
the unpaired nucleon contributes to the vector momentum distribution. On the contrary, 
the scalar momentum distribution receives contribution from the core nucleons as well as 
from the unpaired nucleon, but the latter depends more strongly on the specific structure 
of the particular nucleus under consideration. In the contribution from the core, the A = 
multipole (occurring also for unpolarized nuclei) dominates, while in the contribution from 
the odd nucleon the multipoles with A > are comparable to the A = multipole. Thus, 
it is interesting to focus on transitions in which the knock-out nucleon is the unpaired 
one, and the intrinsic structure of the residual nucleus is basically given by that of the 
even-even core of the parent nucleus. This is the case for transitions to the low lying 
states in the residual nucleus which are populated when the odd nucleon is knocked out. 

We consider a transition from a 100% polarized nucleus in the direction P* = (Q*) to 
a discrete state Jr in the residual nucleus of unobserved polarization. The spin dependent 
transition density in momentum space can then be written as 






J J x^ 

-J J 0; 



E (-l)'"''f_M M' ^]{JM\af,JjRMn){JM'\a+MRMRr , (48) 

M'Mr \ '^7 



X 

MM-ivin 



where we have replaced Eq. (|32D into Eq. (^) and used composition of rotation matrices. 
In addition here we specialize to states Jr in the ground state rotational band {K^ = 0+) 
of the residual nucleus. The matrix elements entering into Eq. (^8]) are in this case 



{JM\a^JjRMR) = ^^Y. j d^T^iflM-^'J^^^) 

X {Vi,^{uj){^,\a^^\%) + {-lY-''Vi_j,{uj){^-,\a^^\%)} (49) 

with <l>o the intrinsic ground state of the even-even core. 

Using the Ij-wave expansion of the single-particle wave functions in the intrinsic frame 
(see Eqs. (||) and (0)), and transforming to the laboratory frame we find that 

($,|a+^|<l>o) = {xj:\Mp,s)r = Y.T.4>i\p){£K-J:l S|j K)Yil{Q,)V%K-Au^) . 

£j me ^ 

(50) 
A similar expression is obtained for (^^lap^l^o) replacing K by —K and multiplying 
by a factor (-1)^"^ in Eq. (|0]). 
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Upon substitution of Eq. ( |50D into Eq. (^) and integration over the direction of the 
intrinsic frame we get 



V^ Pi mfH \ I 



The factor (1 + (— 1)"^^) appears from the sum of the two contributions in Eq. (^Qf), and 
keeps track of the well known property that O"*" bands only contain states with even J 
values, as is the case for the ground state bands of even-even axially symmetric nuclei 
with reflexion symmetry through a plane perpendicular to the symmetry axis. 

Defining scalar and vector transition densities in analogy to the scalar and vector total 
momentum distributions, 

j^iJ^Jn) ^^^ ^*^=rpr (m(^^^«) (p, n*)) (52) 

Mi'-'^^p, Q*) = Tr (m^'-'^\p, fi*)a,) , (53) 

we readily find after substitution of Eq. ( |51| ) into Eq. (^) and computation of the traces 
in spin space: 



A=even \ / 

/ J A 

£ £' 1/2 

Jr J 3 
K -Kj \ K -Kj \J J Jr 



ij £'j' 



(54) 



Mi'-'-\p,Q*) = v^ E E T.i-^ry,\^*)^'-,(^p) 



L=even A=L±1 A* 

X (' ' '](-iy-i' i ^)in 

\a jj, — a —jjij y—J J \}j 

X EEegFwl^' f 1/2I (-i)^'+^'+^(-i)^-^' 

ij er {£ j I/2J 



K -K \ K -K \j J J, 



R\ 



'r 



(Jr J J \ (Jr J y \ \j' J M J2 ^1 + (_i)..) . 

(55) 

Note that these scalar and vector transition densities are made up of the same building 
blocks (the odd nucleon densities Q^'j^ "' (p) defined in Eq.( ^8])) as the intrinsic and labora- 
tory scalar and vector momentum distributions, respectively, and have similar expressions. 
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The only difference is that for each A multipole the sums over ij, i'j' are now restricted 
to angular momentum values satisfying \J — Jr\ < j, j' < J + Jr; \j — Jr\ < f < j + Jr, 
while in Eqs. (P^)-(P7|) and (^Bp-(^) the sums over £j, i'j' do not have these restrictions. 
Furthermore, using the relation 



E J^^' (i 



Jr 

\2.J+X 




f 



-K 




= i-iy-' [{-ly'^^j K j' -K\XO){JKJ -K\XO) 

X 5A,odd(-l)-^+^"+^^(j K f K\X 2K){J K J K\\ 2K)\ (56) 



that we have derived from the relations in Appendix II of Ref. p3[, one can show that 



^M('^'«)(p,fi*)= Y. P^{cose;)nf\p)G{X-,JK), (57) 

Jji A=cvcn 

which is the odd nucleon contribution to the total scalar momentum distribution M (p, Vt*) 
(see Eq. (|36|)), and that 

^ Mi'-'-\p, n*) = Mi'\p, Q*) (58) 

Jr 



We would like to remark that in Eqs. (|57| ) and ( pq) the sum over Jr runs over states 
belonging to the ground state band in the residual nucleus. Therefore, these equations 
tell that measuring the transition densities to each state in the ground state band one can 
map out entirely the intrinsic momentum distribution spin matrix of the odd nucleon. In 
addition, Eq. ( ^8]) tells us that all possible information contained in the vector momentum 
distribution of the polarized target nucleus can be obtained by measuring the transition 
densities to the ground state band in the residual nucleus. The same is true for the odd 
nucleon contribution to the scalar momentum distribution (see Eq. (0)), which as noted 
above is more interesting than the core contribution in the sense that it depends more 
on the nuclear structure and on the polarization. Obviously, transitions to higher excited 
states will bring information on momentum distributions in two quasiparticle excitations, 
vibrational excitations,..., but will not add additional information on the spin dependent 
momentum distribution of the polarized nuclear target. 

As already mentioned, the only difference between the partial and total momentum dis- 
tributions in laboratory frame is the restriction imposed by the Jr value. This restriction 
however, may cause the transition densities to have quite different structures depending 
on the Jr value and on the particular nucleus considered. This is shown in Figures 5 and 
6 where we represent M^^^^^'^p, (]*), M^-^^-^^^p, fi* = 0), and M^-^^^^^ (p, fi* = 0) for 
Jr = and Jr = 2, respectively. The longitudinal (M; ~* ) and transverse (M^ ^ ) 



components (with respect to the polarization direction) of the vector transition densities 
are defined in analogy to M^ and MJ; in Eqs. (^)-(^3|). Fig. 5 is for the transition 

from the ground state in ^^Ne to the ground state in ^°Ne. The scalar transition density 
is shown on top, and the longitudinal and transverse components of the vector transition 
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density are shown below. For this transition {J -^ Jr = 0) the momentum dependence 
of the transition densities does not change appreciably in going from ^^Ne to ^^Ar, only 
their strength change. Thus, Fig. 5 also serves to represent the transition densities for the 

ground state in ^^Ar to the ground state in ^^Ar. For this purpose one has to multiply 
the results shown in Fig. 5 by a factor ^ 14, which is the ratio between the ^3/2 strength 

in the odd nucleon wave functions of ^^Ar and ^^Ne. ^ ^ 

On the contrary, the momentum dependence changes much in going from ^^Ne to ^^Ar 
when considering the transition to the first 2+ states of their respective residual nuclei 
(^''Ne and ^^Ar). This is clearly seen in Fig. 6, where the plots on the left correspond 

to 2iNe (3/2+) ^20 j^g (^2^) and the plots on the right correspond to ^'^Ar (3/2+) ^^^ 
Ar(2]^). Scalar momentum distributions are shown on the top and longitudinal (/) and 
transverse (t) components of the vector momentum distributions for these transitions are 
shown below. As already noted when comparing the intrinsic momentum distribution for 
^^Ne and ^^Ar in section II, in the case of ^^Ar the odd nucleon wave function contains 
mainly a mixture of s and d waves, whereas in the case of ^^Ne it contains mainly d 
wave. For the transitions J ^ Jr = only the ^3/2 wave component contributes in either 
case, thus resulting in similar momentum dependence of the transition densities. For the 
transitions J — ;► Jr = 2 all positive parity waves with j = ^, |, f , | may contribute, but 
the j = I is only present in the odd nucleon of ^'^Ar, not in ^^Ne. This, together with the 
fact that also the |, |, | components enter with different amplitudes in each case, explains 

why the momentum dependences in Fig. 6 are different for ^^Ne and '^^Ar. Finally we 
would like to remark that the transition densities shown in Figs. 5 and 6 are the dominant 

contributions to the total scalar and vector momentum distributions in ^^Ne and '^^Ar. 
Transitions densities to higher excited states 4+, 6+,... in the ground state band of the 
residual nuclei are much smaller than the ones seen in Figs. 5 and 6. 

III. 3 Comparison with Spherical Case 

It is interesting to compare the results in previous sections with the ones corresponding to 
spherical nuclei. The case of spherical nuclei can be thought of as a limiting case when the 
deformation of the mean field goes to zero. In this limit the single particle wave functions 
are eigenstates of angular momentum and the system does not have a preferred direction 
unless the nucleus is polarized by an external field. Hence the main two differences with 
the deformed case are that now, i) each single particle wave function has a single £j-wave 
component {(f>lj{p) = 61^1-6 jj.Rn^i.j-{p)) , and ii) there is no distinction between intrinsic 
and laboratory frames. We denote by {i} = {riiiiji} the set of orbitals filled by the 
nucleons of the even-even core and by {k} = {rikikjk} the odd nucleon orbital, which has 
jk = J and ik fixed by the parity and J values. 

The scalar momentum distribution for the fully polarized nucleus is then given by 

M^'^''"'-ip,n^)= Y: PA(cos^;)(5,,on^P'^-(p)+n°/'^'^P^-(p)) (59) 

A=cvcn 
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where "sph." stands for spherical hmit; the odd nucleon multipoles are given by 



n 



odd, sph. 
A 



ip) 



-k^ Jk I 6 



An 



ik ik A^ 

0; 



Rn^MlTn n nl(-l)^'="- 



(jk jk M ijk jk A| 



{Jk 



-Jk 0/ 14 4 M '"^ 



and n^'P^' {p) is the spherically symmetric momentum distribution of the core 



n 



^P^-(p) = no(p) = ^E%'(2J. + mn.iaAp)\' ■ 



(60) 



(61) 



A thorough study of the properties of n{p) for even-even spherical nuclei can be found 
in Refs. [^, |2^ . The differences between scalar momentum distributions in the spherical 
and deformed case are now made explicit by comparing Eqs. (p9D-(|61|) with Eq. (p6D (see 
also Eq. (0)). The core has now only a monopole contribution, while in the deformed 
case the core has also small multipoles n\{p) with A > 0, that approximately average to 
zero after integration on p (see Fig. 2 and table II). 

For the odd nucleon contribution the sum over ij, i'j' in Eq. ( PTj ) is now reduced to 
a single term ij = i'j' = ikjk, with i^jk the values corresponding to the spin and parity 
of the nucleus, and with Vtk = K = jk = J- Therefore, in the spherical case all the 
A-multipoles have the same p dependence, whereas in the deformed case for each A value 
we may have a different dependence on p. The geometrical coefficient G{\\ JK), which 
appears in the deformed case from the transformation between intrinsic and laboratory 
frames, is not present in the spherical case. This reflects the fact that now, "a priori", 
there is no internal preferred direction. 

The longitudinal (/) and transverse (t) components of the vector momentum distribu- 
tions are 



?(J)sph. 



I * I 



x(- 



0) 



\R 



riktkJk 



Xp)P 



47r 



jk,j 2^ z^ 

L=even A=L±1 




(L A 1] 




h jk \ 


• • 


[4 Jk \\ 






p^cose; 



(62) 



where for simplicity we have chosen the polarization direction along the laboratory z axis. 
Again, by comparing Eq. (p^ ) to Eqs. (^2|) and (P3| ) (see also Eqs. (^) and (PH)) we 
see that the differences with the deformed case are the absence of the geometrical factor 
G{\; JK) and the restriction of the sum over ij, i'j' to a single term with i^jk dictated by 
the spin and parity of the nucleus, and with Qk = jk = J- It is important to remark that 
this restriction produces equal p-dependence of the various scalar and vector momentum 
multipoles. 

It is also interesting to remark the similarity between Eq. (|6^) , that gives the longi- 
tudinal and transverse vector momentum distributions with respect to the polarization 
direction, and Eqs. ( ^5]) and (|26|) , that give the intrinsic vector momentum distribution 
of the deformed nucleus along the symmetry axis and in the perpendicular plane. If one 
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restricts Eq. ( pTf ) to a single term with ij = i'j' = ikjk, and takes Qk = jk = J, the above 



mentioned expressions are identical to Eq. (62). This reflects the fact that for spherical 



nuclei polarizing the nucleus is the only way to define an internal preferred direction for 
the system. 

Expressions for transition densities in this case can also be easily derived. Obviously in 
this case to map out the vector momentum distribution and the odd nucleon contribution 
to the scalar momentum distribution one needs only to consider the transition in which 
the residual nucleus is left in its ground state (J/j = 0). It is easy to show that in this 
case the vector momentum distribution for the transition J ^ J^ = is also given by 
Eq. (0), and the scalar momentum distribution is given by the odd nucleon contribution 
in Eq. (^). For the case of a spherical nucleus with odd nucleon in a single rfs orbital 
these momentum distributions have the same structure as the one shown in the top plot 
of Fig. 4. 

IV. Summary and Final Remarks 

We have studied spin dependent momentum distributions in deformed nuclei on the basis 
of the selfconsistent mean field approximation. The spin degree of freedom is relevant 
for odd-A nuclei that can be polarized. We consider the case of axially and reflection 
symmetric mean fields and discuss momentum distributions in the nuclear ground state. 
We first study spin dependent momentum distributions in the intrinsic frame for a set 
of deformed nuclei in the sd-shell (^^Ne, ^^Ar, ^^Mg), and then relate them to momen- 
tum distributions in the laboratory frame and to transition densities in momentum space 
that can in principle be measured in one nucleon knock-out reactions. We do this in 
a systematic way by decomposing the spin dependent two by two matrix M^^i{p) into 
scalar and vector components in spin space. We show that the scalar and vector momen- 
tum distributions that can be measured in the laboratory are intimately related to the 
corresponding scalar and vector intrinsic momentum distributions. Thus, measuring the 
former gives information on the latter. 

The study carried out for the selected nuclei (^^Ne, ^^Ar, ^^Mg) shows that the even- 
even cores, which contribute only to the scalar momentum distribution, are dominated 
by the A = multipole but contain also quadrupole (and higher) multipoles that are 
small and tend to average to zero upon integration on p. Actually, this isotropy condition 
((pi) = (Pz)) that was already discussed for the case of even-even nuclei in Refs. [^ |^, 
is found here to be also approximately satisfied for the odd-A deformed nuclei. More 
involved structure is found in the vector momentum distribution and, to a lesser extent, 
in the odd nucleon contribution to the scalar momentum distribution. 

The vector momentum distributions are found to have very rich structures, particu- 
larly for ^'''Ar which is the nucleus with lowest K value in its ground state. The richness 
in structure of the vector momentum distribution, as well as of the odd nucleon contribu- 
tion to the scalar momentum distribution, is found to decrease with increasing K values. 
Information on this internal structure can be gained by measuring the transition densities 
in momentum space by one nucleon knock-out reactions from a polarized target. It is 
also shown that most of the interesting information can be obtained from transitions to 



the lowest states (ground state band) in the residual nucleus. An application to A(e, e'n)B 
reaction on ^^Ne has already been madeEJ and applications to the same type of reaction 
on ^^Ar and ^^K are now under studyO. The expressions relating intrinsic and labo- 
ratory momentum distributions have been derived here within a similar philosophy to 
that in Refs. |jT6|, ^ concerning the relations between intrinsic and laboratory multipoles 
measured by (e,e') reactions. 

We have also compared the results on deformed nuclei with results for spherical nuclei. 
This comparison shows that in the spherical limit the structure of vector and scalar mo- 
mentum distributions is much simpler, even when we consider a polarized odd-A nucleus. 
Obviously for unpolarized nuclei, whether spherical or deformed, only the monopole part 
of the scalar momentum distribution enters into the picture; neither the vector momentum 
distribution nor the higher multipoles of the scalar momentum distributions are accessible 
in unpolarized nuclei. 

The results presented have been obtained with the SKA interaction, similar results 
are obtained with other Skyrme type forces. It would be interesting to check whether the 
results presented here are modified when finite range effective interactions, like Gogny 
forceEj, are used. In principle this would allow to study effects of short-range dynamical 
correlations not included in the present calculations. It would also be interesting to study 
the results obtained with relativistic deformed HF calculations. The latter calculations 
may be superior to the non-relativistic ones for the study of spin degrees of freedom. 
In particular for odd-A nuclei the relativistic formalism allows to take into account in a 
simpler way effects due to time reversal non invariant terms in the mean field that have 
not been considered in the present work. In the future we plan to study these points. 
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Figure Captions 

Figure 1: Contributions from tlie even-even core (left) and from tlie odd nucleon (riglit) to 
tlie scalar momentum distribution in ^^Ne (top) and ^^Ar (bottom). 

Figure 2: Monopole (no) and quadrupole (n2) contributions (in fm^) to the even-even core 
momentum distributions (normalized to 1) of ^^Ne and '^^Ar from DDHF calculations 
and from Nilsson model calculations. The results obtained with the Nilsson model 
without major shell admixtures are labelled with primes. 

Figure 3: Components of the vector momentum distribution parallel and perpendicular to the 

symmetry axis for ^^Ne (top) and '^'^Ar (bottom). 

Figure 4: Core (bottom) and odd nucleon (top) momentum distributions in ^^Mg. The vector 
momentum distribution has a single component parallel to the symmetry axis that 
coincide with the odd nucleon contribution to the scalar momentum distribution 

(M,(p)=n°<i'i(p)). 

Figure 5: Scalar (M ^ {p,il*)) and vector components {Mij^'^^\p,VL* = 0)) of the tran- 
sitions densities in momentum space (in fm^) for J ^ Jr = in ^^Ne. The results 
for ^^Ar are similar except for a scale factor ~ 14. 

Figure 6: Same as Fig. 5 for the transition J -^ Jj:j = 2 in ^^Ne (left) and in ^'^Ar (right). 
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Table I: Results of DDHF calculations for binding energies, proton and mass quadrupole 
moments, r.m.s. radii, moments of inertia and decoupling parameters. 



Nucleus B (McY) Q"" {iia^) Q^ (fm^) (r^)2 (fm) J (MeV-i) 



2iNe 


-165.1 


44.0 


92.1 


3.00 


4.19 


20Ne 


-155.5 


44.0 


86.8 


3.00 


4.38 


37Ar 


-312.3 


-43.1 


-74.6 


3.41 


1.89 


36 Ar 


-301.9 


-43.1 


-84.8 


3.41 


2.50 



-1.159 
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Table II: Values of quadrupole deformation parameters in p-space and in r-space from Nilsson 
model calculations with (AA^ 7^ 0), and without (AA^ = 0) major shell admixtures, 
compared to results of DDHF calculations. 





2iNe 


20Ne 


37 Ar 


36 Ar 




(3P /?'■ 


pP p'' 


pP /9^ 


pP /?'■ 


DDHF 


0.041 0.402 


0.050 0.405 


0.014 -0.145 


-0.011 -0.170 


Nilsson^^^o 


0.022 0.445 


-0.002 0.487 


0.028 -0.138 


0.012 -0.171 


Nilsson^^=o 


0.232 0.232 


0.234 0.234 


-0.062 -0.062 


-0.084 -0.084 
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